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Abstract. Decomposing the domain of a function into parts has many 
uses in mathematics. A domain may naturahy be a union of pieces, a 
function may be defined by cases, or different boundary conditions may 
hold on different regions. For any particular problem the domain can 
be given explicitly, but when dealing with a family of problems given in 
terms of symbolic parameters, matters become more difficult. This article 
shows how hybrid sets, that is multisets allowing negative multiplicity, 
may be used to express symbolic domain decompositions in an efficient, 
elegant and uniform way, simplifying both computation and reasoning. 
We apply this theory to the arithmetic of piecewise functions and sym- 
bolic matrices and show how certain operations may be reduced from 
exponential to linear complexity. 



1 Introduction 

The goal of this paper is to develop general methods to work with domains 
having symbolically defined parts. The raison d'etre of symbolic mathematical 
computation is to compute and reason about general expressions, rather than 
working with particular values valid only at specific points. Matters are simplest 
when variables range over a domain of interest and all expressions are valid over 
the entire domain. Sometimes it is useful to perform simplifications or other 
operations that arc valid over part, but not all, of the domain. In this situation, 
software systems may or may not record the excluded region. But we are not 
always so fortunate to have this one-region situation. More generally, the domain 
of interest may be made up of several pieces with expressions taking different 
forms on different parts. Moreover, the demarcation of the parts may be defined 
symbolically. This paper explores how to represent such expressions concisely in 
a uniform way that simplifies computation and reasoning. 

When we do arithmetic with piecewise functions defined on explicit parti- 
tions, we can do a mutual refinement of domain partitions to obtain regions that 
may each be handled uniformly. When the partitions are defined symbolically, 
however, we obtain a massive explosion of cases — for N binary operations on 
functions of k pieces there are fc^ potential regions. We say "potential" regions 



because, of this large number, not all of the regions are in fact feasible. Fur- 
thermore, it is usually not possible to determine which regions are feasible and 
which are not. For example, the sum J2i^=i fi{^) of functions 




for a; < ki, 
Ai for X > ki. 



has X]i=o-^'/^' possible orderings of the fc^, with each ordering having between 
2 and A'^ + 1 regions. There is no ordering in which all 2^ regions are realized. 

We take the view that it is generally preferable to have a single compact 
closed-form expression rather than a collection of cases, even if it means intro- 
ducing some new operations. For example, we are perfectly satisfied using the 
Heavisidc step function and giving the sum of the as 'Y^=\ AiH{x — ki). 

In this paper we show how hybrid sets, a variation on multisets allowing 
negative multiplicities, enable us to write elegant closed form expressions of the 
form we desire. This use of hybrid sets also allows us to define a generalised notion 
of partition, where symbolically defined parts are combined in more useful ways 
than the usual set operations. Our approach unifies and generalises a number of 
other techniques, such as the use of oriented regions for domains of integration. 

Introducing new operators or generalising existing ones to write single closed 
form expressions is more than just a cosmetic re- arrangement. It allows one 
to perform arithmetic and simplifications on whole expressions, and to reason 
about the expression and about the regions themselves. It is already customary 
to do this for certain operators. For example, by defining fdx = — fdx, 
it becomes possible to write identities that hold universally. Then we have that, 
independent of the relative order of a, b and c, and subject to / being defined 
on the requisite domains, 

f fdx= f fdx+ f fdx. (1) 

J a J a J c 

With a little work, we can also generalise the integral formula to integrating 
over oriented subsets of M". Some authors similarly adjust the definitions of 
other operators to obtain universally true statements. Similarly, Karr [8] defines 
the summation operator Em^i<„ so that Em^i<„ = - T,n^i<m ^^^en n < m. 
This allows equations such as the following to hold for any ordering of £, to, n: 

^ (g{^ + l)-gi^))=g(n)-gim), ^ /(») = ^ /(*) + 

This paper formalises and extends these ideas in several ways, giving a gen- 
eralised framework for domain partitions and piecewise defined functions. We 
first introduce some preliminaries and hybrid sets (§2) and then their general- 
isations to our notions of generalised partitions and hybrid functions (§3). We 
then present how we can decompose domains of hybrid functions to allow for 
the combination of their symbolically defined pieces (§4), before presenting some 
applications (§5) and discussing some concrete examples (§6). 



2 Preliminaries 



2.1 Partitions and piecewise functions 

The domain of definition of various matliematical objects (functions, vectors, 
matrices, sequences, etc) may naturally be decomposed into a (disjoint) union of 
pieces where our object is then defined uniformly. When the pieces are given as an 
explicit union of provably disjoint sets which form a partition of the domain, the 
interpretation of a given expression is reasonably straightforward. More formally, 

Notation 1 We use the notation Ci^jXi, or, more briefly, CiXi to describe a 
collection of elements Xi,X2, ■ ■ ■ , indexed by a set I. For n e N, we denote by 
[n] the set {1, . . . , n}. 

Definition 1. A partition of a set U is a collection CjPi of pairwise disjoint 
sets such that [Jj Pi = U. 

A partition, CjPi, induces a total fimction X : U ^ I which gives the index of 
P* where each u ^ U sits. Piecewise expressions are then defined on top of a 
partition. 

Definition 2. A piecewise expression over a set U is a collection Ci{Pi,ei) 
where CjPi is a partition of U and each Ci is an expression. 

Typically, each Cj contains a distinguished variable y which is interpreted to 
range over U. 

Definition 3. We say that f : U ~> S is a picccwise-defined function if we have 
a collection Ci{Pi, fi) where CjPi is a partition of U , "ii € I. fi : Pi ^ S, and 

^x:U.f{x) = fx(:,){x). 

We call Ci{Pi, fi) the definition of f, and each fi a piece of f. 

It is important to note that piecewise-defined functions use their argument in 
two different ways: once geometrically by choosing a set Pi "over" which to 
work, and once analytically to evaluate a function /j. The definitions above are 
straightforward generalisations of those found in [5]. 

Lastly, we have that arithmetic operations on piecewise-defined functions 
operate component- wise. Note that we will silently use the convention that op- 
erations defined on the codomain of a function are lifted pointwise to apply to 
functions, in other words (/ + g){x) = f{x) + g{x). 

Proposition 1. Let f,g : U S be two piecewise-defined functions on the 

same partition C I Pi of U , withCjfi (respectively Cjgi) the collections of pieces 
of f (resp. g). Further, let -k : S x S ^ S. Then Cj {fi * gi) is the collection of 
pieces of f -k g over the partition CiPi . 



Note how the partition is entirely untouched. This "separation of concerns" is 
what enables us to separate the issues of domain decompositions from arithmetic 
issues of piecewise-defined functions (and expressions). 

To simplify our presentation, we introduce a domain restriction operation and 
a join combinator on (partial) functions, to allow us a more syntactic method 
of "building up" piecewise functions. These are quite similar to Kahl's table 
composition combinators [7]. 



Definition 4. The restriction of a function f to a domain specified by a set 
A is 

' f{x) ifxGA 
_L otherwise 



f\x) ::= 



Definition 5. The join, f (D g, of two (partial) functions f and g, is defined as 
{f(Dg){x) ::= < 



f{x) if f(x) is defined and g(x) is undefined 
g{x) if g(x) is defined and f(x) is undefined 
_L otherwise 



This allows us to rewrite a piecewise-defined function / defined by Ci{Pi, fi), in 
terms of its pieces as 

f^f['G)f['Q...Qf^-. 

But oiu- goal is to work with piecewise-defined functions where we have a 
symbolic partition. We need some new tools for this, which we will develop in 
the next two sections. 



2.2 Hybrid Sets 

We consider an extension of multisets, in which elements can occur multiple 
times, to hybrid sets, where the multiplicity of an element in a hybrid set can 
range over all of Z, instead of just No- Thus a hybrid set, over an underlying 
set U, is a mapping U ^ Z, i.e., it is an element of . We use the following 
definition, adapted from [9]: 

Definition 6. Given a universe U, any function H : U Z is called a hybrid 
set. 

We can immediately define some useful vocabulary for working with hybrid 
sets. 

Definition 7. The value of H{x) is said to be the multiplicity of the element 
X. If H{x) ^ 0, we say that x is a member of H and write x G H; otherwise, 
we write x ^ H. The support of a hybrid set is the (non-hybrid) subset S of 
U where s € S s (E H; we will denote the support of H by suppiJ. We 

(re)use to (also) denote the empty hybrid set, i.e. the hybrid set for whom all 
elements have multiplicity 0. 



Notation 2 We use the notation {\x^^ ,x^^, . . .|} to describe the hybrid set con- 
taining elements Xi with multiplicity mi, X2 with multiplicity m2, etc. While our 
notation allows writing hybrid sets with multiple copies of the same element with 
different multiplicities, these denote the same hybrid set as that denoted by the 
normalised form with one copy of each element with a multiplicity which is the 
sum of the multiplicities of the copies of that element in the non-normalised 
form. Thus {\a^ , b\ , b'^\} = {\a'\b^\}. 

Set unions are usually defined by the boolean algebra structure (and more 
specifically, via V) of the membership relation. For hybrid set, this is replaced 
by arithmetic over Z. 

Definition 8. We define the sum, A(B B of two hybrid sets A and B over a 
universe U, to be their pointwise sum. That is {A®B){x) = A{x) -\-B{x) for all 
X G U. We similarly define their difference, AqB to be their pointwise difference, 
and their product, A ® B to be their pointwise product. Let QB denote _B. 

In other words, we do not use operations AUB, AUB and A\B for hybrid sets, 
but just ®, e and 0. We can easily establish some identities such as {A(BB)qA = 
B, AeA = <!) and A ® (qB) = AqB as these follow directly from Z. 
Putting all of this together, wc get: 

Proposition 2. Z''^ is a Z-m,odule. 

Proof. The abelian group structure is given by (Z^,©,0,0), and Z acts on 
hybrid sets by nH = n - H{u). 

We need two more technical definitions, which will be useful later. 

Definition 9. We say that two hybrid sets A and B are disjoint if A^ B = $. 

Definition 10. We call a hybrid set reducible if all its members have m,ultiplic- 
ity 1. We define a reduction function, 7?,(-), on reducible hybrid sets that returns 
the (normal) set of members of the hybrid set. 

We should note that these hybrid sets (sometimes also called generalised sets) 
have been studied before. Hailperin [6] makes the case that Boole [3] actually 
started from hybrid sets for his algebraization of logic, but restricted himself 
to nilpotent solutions of the resulting equations, which then correspond closely 
to our modern notion of Boolean algebra. Whitney wrote two nice papers [13, 
14] taking up the theme of algebraising logic via characteristic functions. He 
does allow arbitrary multiplicities, and derives some nice normal forms for cer- 
tain kinds of partitions, in a way foreshadowing some of our own results (see 
§3.2 and §4). Blizard [1] focuses on multisets (disallowing negative multiplici- 
ties) but has an extensive bibliography of related works, several of which being 
on (mechanised) theorem proving; he then formalised sets with negative mem- 
bership in [2] . Blizard concentrates on concepts of union and intersection which 
closely resemble those of normal set theory, although he does also define the sum 
union (but not other related concepts). Burgin [4] lists several more works on 
hybrid sets, some reaching back to the early middle ages. Syropoulos [12] gives 
a very readable introduction to both multisets and hybrid sets. 



3 Generalisations 



We now revisit a few basic mathematical constructs and show how they may 
be modified to work with hybrid sets. This will provide the machinery that we 
need for symbolic domain decomposition. First, we will examine the notion of a 
hybrid function on a domain. We then show how sets and hybrid sets may be 
decomposed using a notion of generalised partitions — an extension of parti- 
tions to the hybrid set case. We then address the practical issue of how to make 
two hybrid partitions compatible by constructing a common refinement. Finally, 
when working with functions defined over hybrid partitions, we need some way 
to compute values. Over any given point in the domain, we need to know which 
functions must be evaluated in computing the final value, which is rather com- 
plex for hybrid functions over generalised partitions. For this task, we introduce 
the notion of pseudo-functions. When expressions on generalised partitions are 
evahiated, these pseudo-functions avoid evaluating at places where the functions 
are undefined or where the values are not needed. This allows us to deal with the 
cases, as in equation (1), where component functions are not defined on some 
parts of the domain decomposition but any application of the function in those 
places would anyway have multiplicity zero {i.e. not be used). 

3.1 Functions of hybrid domain 

It turns out that a useful definition of a "function" involving hybrid sets is not 
entirely straightforward. Defining its graph is easiest. The underlying intuition 
is that we capture the restriction of a function to a domain through the mul- 
tiplicities of the elements of the function graph in a hybrid set. The hybrid set 
of a single element of the function graph for element a; in t/ is of the form 
Therefore the scalar multiplication of that set by the multiplicity 
of a; in A will impose the appropriate restriction. We use our function restriction 
notation of Def. (4) only for this hybrid version of function restriction henceforth. 

Definition 11. Let A he a hybrid set over U , B C U , S a set and f : B ^ S a 
(total-on-B) function. A hybrid function f^:UxS^Zis defined by 

f^ = ^Aix){\{x,f{x)y\} 
xeB 

Note how the hybrid set Z-module structure automatically takes care of re- 
stricting the sum over the support of A. Caution: some hybrid sets do not form 
a hybrid function (for example {|(1, 1)^, (1, 2)^|} is not a hybrid function). 

Our definitions work just as well with partial functions as with total fimctions. 
But if for a hybrid function , / is undefined at some point in the support of 
F, will not be defined at that point either. So, without loss of generality, we 
can always restrict F to where / is defined. For the remainder of this paper, we 
shall assume this, i.e. whenever we write a hybrid function f^, f is total over 
supp F. 



Definition 12. We call a hybrid function reducible if the hybrid set H is 
reducible. We extend TZ{-) in this case by 

' \± ifH{x) = 

We can generalise the join combinator to hybrid sets. This definition is quite 
central to "making things work" . 

Definition 13. The join, f^ ® , of two hybrid functions f^ and (with 
codomain B ), gives a hybrid relation, a subset ofU x B x Z given by 

f^Og"" ::= reff^ 

This is a rather "dangerous" definition, as it moves us from the land of functions 
to that of relations. In other words, it is quite possible that f^ (I)g'^ restricted to 
[/ X B is no longer the graph of a function, but the graph of a relation. But this 
extra generality will be quite useful for us, although we will have to prove that in 
the cases which interest us, the resulting hybrid relations are in fact (reducible) 
hybrid functions. 

Theorem 1. Let A, B be hybrid sets over U, S an arbitrary set, and f : U ^ S 
a total function. Then 

1. TZ{f^) is the empty function, 

2. f^'(J) fA = /2'4 

3- f"^ CD f^ = f^®^, and thus a hybrid function, 

4. For g : U ^ S another total function, then f^ d) g^ ^ {f d) g)^®^ if and 
only if A® B — % (where f d) g is the join of regular functions) . 

5. Let Hi , H2 be hybrid sets, with supp Hi and supp H2 disjoint, fi : supp Hi — >■ 
S and f2 : suppiJa ^ S, then /f ^ (D /f ^ = (/i d) /a)^^®^^ 

The proofs are omitted, and follow straightforwardly from the definitions. Note 
the strong dichotomy between (3) and (4), which comes from the fact that the 
non-hybrid ® is designed to work with functions defined over separate regions. 



3.2 Generalised Partitions 

Theorem 1 tells us that some collections of hybrid sets are better than others. 
Being disjoint is much too strong a property. Nicely, for hybrid sets, partitions 
easily generalise in useful ways. 

Definition 14. We define a generalised partition of a (hybrid) set, P , to be a 
finite collection of hybrid sets, C[n]Pi, such that Pi ® P2 © • • • ® -Pn = -P 

All set partitions of a set are also generalised partitions. Conversely, a generalised 
partition of a reducible set is a set partition if and only if each generalised 
partition element is reducible. 



Remark 1. We have lifted the disjointness condition on partitions. For some- 
thing to be called a partition of P, it is necessary that the result be equal to P. 
Still, P belongs to a larger universe U, and a generalised partition's pieces range 
over U. As long as, in the end, all elements of U\P have multiplicity 0, we get 
a generalised partition. In this way, we have also lifted the coverage condition. 

Proposition 3. For any generalised partition C^n]Pi of a hybrid set P over U, 
arbitrary set S, and any function f : supp P ^ S, 

= /^i ® . . . ® fPn = ;PieP2e...eP„ 

is a hybrid function. 

For brevity, we sometimes write f^ for either of the right-hand side expressions 
above. When we want to join different functions over a partition and still get a 
function, wc have to be careful and ensure we arc joining "compatible" functions. 

Definition 15. Let Pi,P2 be a generalised partition of a hybrid set P over U, 
S an arbitrary set and f^^ : Pi — )• 5, g^^ : P2 ^ S hybrid functions. We 
say that Pi,P2 is a compatible partition for f,g if f{x) = g{x) for all x G 
supp Pi n supp P2 ■ 

Theorem 2. Using the same notation as above, f^^ (D g^^ is a hybrid function 
if and only if Pi, P2 is a compatible partition for f, g. 

It is important to note that although ® is an associative, commutative oper- 
ation, the notion of compatibility, while commutative, does not lift to a simple 
associative condition. 

Remark 2. Some of our computations will purposefully use incompatible parti- 
tions. Note that 

(/^ CD 5^) (D 5®^ = f(D {9" CD fl®^) = f(D9' = f 

but that ® g^ is in general a hybrid relation, yet the final result is a hybrid 
function whenever is. We will "design" our hybrid partitions with this feature 
in mind. 

3.3 Refinement 

To do arithmetic with hybrid functions, we first need the notion of a refinement 
and a common refinement. This is similar to the treatment of [5] for piecewise 
functions. 

Definition 16. A refinement of a generalised partition CjPi of P is another 
generalised partition CjQj (of another hybrid set Q not necessarily equal to P) 
such that for every i G I there exists a sub-collection Qj^ of Qj such that Pi = 
Qji ® Qj2 . . . ® Qj^ ■ A common refinement of a set of generalised partitions 
is a generalised partition that is simultaneously a refinement to every partition 
in the set. 



A refinement in this sense may well seem "larger" than the original partition, as 
in the next example. 

Example 1. Let the interval P = [0,1], seen as {l^*^!}, and Ii = [—1,0), I2 = 
(1,2], I3 = [-1,2], then Q = {|7f \ 72~S -^3 1} is a refinement of P. 

Definition 17. A refinement is called strict if, for each generalised partition 
being refined, the support of the associated sub- collection is equal to the support 

of the generalised partition it refines. 

Example 2. {|[0, 1]^|}, {|(1,2]^|}, ||(2, 3]^|} is a common strict refinement of the 
two (trivial) hybrid partitions {|[0,2]^|} and {|(1,3]^|}. 

3.4 Pseudo-functions and pseudo-relations 

As seen in the example above, a refinement may "spill over" the original domain, 
so that if we look at a hybrid function f^ where the underlying / is defined 
exactly on (the support of) P, f'^ evaluated "pointwise" will not make sense. 
Nevertheless, wc want f^ = f^ ■ To achieve this, we apply the lambda- lifting 
trick already used in [5]. 

Definition 18. Using the same notation as in Def. 11, we define a pseudo- 
function as 

f^^^A{x){\{x,fY\] 
xeB 

i.e. as a member of U x (/7 — > 5) — >■ Z. ^ pseudo-relation is defined similarly. 
The evaluation of a pseudo-function (resp. relation) is defined by mapping each 
point {x,fY to {x,f{x))^. 

The usefulness of pseudo-functions comes from the following property. 

Proposition 4. For all refinements Q of the generalised partition P, f^ — f^. 

In other words, even though the pieces of Q might "spill over" , if we first "sim- 
plify" /"^ by performing 0^ Qi to get P, we get a result f^ which can then be 
safely evaluated. We will elide the ~ to lighten the notation whenever this would 
not lead to confusion. 

Another useful property of pseudo-functions is that in some cases we can 
simplify them, regardless of what the underlying function does. 

Proposition 5. f^ (Dg'^ (D g^'^ = f^ 

One of the chief advantages of pseudo-functions is that we can do some symbolic 
manipulations of expressions in terms of these functions as if they were defined 
on a much larger domain, as long as the eventual term we evaluate does not 
involve any of these "virtual" terms. 

To aid in such computations, when we have pseudo-functions f^ and , 
with /, g : supp P — )• 5, and some binary operation -k : S x S ^ S, we will allow 



ourselves to write expressions such as ★ in the induced term algebra over 
pseudo-functions. As usual, we lift evaluation pointwise, {f*g){x) = f{x)-kg{x). 
Furthermore we say that f = g over a set supp P whenever e supp P.f{x) = 
g{x). In other words, we use extensional equality for the intensional terms / 
and g. This means that properties like commutativity, associativity and having 
inverses lift to the term algebra. As an example, we have that 

Proposition 6. If -k : S x S ^ S is associative and commutative then 



4 Hybrid domain decomposition 

We now have all the ideas necessary to decompose hybrid domains. An elegant 
consequence of the formalism is that it allows us to use linear algebra to construct 
the partitions that we require. 

Let C[n]Ai and C[„L\Bj, be generalised partitions of U. We want to find a 
generalised partition of U that is a common strict refinement of Ai and Bj and 
has minimal cardinality. The cardinality restriction is to minimise the number 
of terms required for a symbolic representation of the resulting domain decom- 
position. Thus we want to choose a minimal generalised partition, C/Pj of U 
such that, in the Z-module of hybrid sets and for some integers aij, bij we have 
®iP^ = U and Vi : l..n. 0^- a,,jP, = A, and Vj : l..m. 0, &»,jP«'= B^. 

Since this forms a system of n -|- m -t- 1 simultaneous equations, of which 
only n + m — 1 can be independent, because A and B are, separately, partitions 
of U, wc need that number of independent variables to solve the system. Thus 
the cardinality of the minimal partition is n -|- m — 1 in the general case, and 
can only be smaller in specific cases if there are some extra dependencies among 

U,Ai,..., An-l,Bi, . . . , Bm-l- 

This result generalises to a decomposition of U into a minimal generalised 
partition that is a common strict refinement of r generalised partitions of car- 
dinality ni,...nr respectively: The minimal partition required, assuming full 
independence of the individual domain partitions, has cardinality 

If all the individual partitions have the same cardinality, n, this reduces to 
r(n- 1) + 1. 

We can automatically compute suitable minimal strict refinement partitions 
U as follows. If we remove the equations for An and B^ from the system of 
equations in order to get an independent set of simultaneous equations, we get 



n + m — 1 equations that can be written as a linear system in the Z-module: 



/ Pi 



C- 



+m-lj 



\Bm-l/ 



( 1 

ai,i 



where C 



I 

ai,2 



^1,1 ^1,2 



1 \ 

^l,n+m— 1 

^n— l,n+m— 1 
&l,n+m-l 



V6 



'm-1,1 Om-1,2 



^m— l,n+m— 1 



Note that C is an integer matrix. Further, the partition choice matrix, C, 
must be invertible and must be an integer matrix so that we obtain integral 
partitions of each domain piece with respect to owx partition P of U . An integer 
matrix is invertible and has an integer matrix inverse if and only if it has a 
determinant of +1 or —1. Hence our problem of constructing an appropriate 
partition reduces to choosing an integer matrix of the form of C such that its 
determinant is ±1. Finally, note that this directly generalises to an arbitrary 
number of piecewise functions, each of an arbitrary number of pieces. 

If we restrict ourselves to triangular matrices, wc can choose C to be any 
integer triangular matrix (upper because the first row of C is all Is) for which 
the product of the diagonal elements is 1. Again, a simple way to do this is to 
choose C to be all Is along the top row. Is along the diagonal and everywhere 
else. Another possibility is all Is in the whole upper triangle. 

For example, two suitable choice matrices with their inverses are 



/I 



1 ... 



Vo 



■•• 

1/ 



/l-l.... 

1 ... 



-1\ 



Vo 



■•• 

1/ 



/I 



/l-l ... \ 
■■• ■•• ■•• : 
: '•• '•. '•. 



Vo... 1/ 



Vo 



•. -1 
1/ 



5 Applications 

5.1 Arithmetic on Piecewise Functions 

We are now ready to generalise the arithmetic properties (see prop. 1) of hybrid 
functions and pseudo-functions. 

Proposition 7. Let Cyn\Pi be a partition of P, = /f^ CD ... CD f^" and = 
Qi^ © . . . ® g^" be two hybrid functions on P ^ S. Let * : {S x S) ^ S, then for 
all X e supp P, 



f{x)^g''{x) = {f^ix)*gi{x)f- ® . . . ® (/„(a;) ★SnW)^"- 



Note how we can apply the above proposition to *g^, by first restricting F 
and G to be over a common support (as x*_L = -L*y = -L), then taking a 
common strict refinement of (the restricted) F and G. 

We would like to lift this strictness condition. We can almost do this with 
pseudo- functions - and with the help of a marked ®, we can. 

Definition 19. We can mark a CD with a binary operation -k : S x S ^ S, 
denoted (D* ■ We define evaluation of (D* by 

{F (S)*f){x) = {F{x) + G{x))^{x,{fk~g){x)f\^ 

This operation clearly inherits properties of -k like commutativity, associativity 
and invertibility. Unlike ®, ®* will always result in a (pseudo) function. 

Proposition 8. Let P* = C[„]-P; be a partition of P, Q» = C^^-^Qj another 
partition of P, and — CkRu a common refinement of P^ and Q*. Let f^ = 
f^^ CD ... (D fn" O'lT-d ~ gf^ (D ■ ■ • CD^m™ be two pseudo functions with codomain 
S. Let -k : (S X S) ^ S , be associative and commutative, then ★ distributes over 
the partition i?, in terms of ®* . By the results of section 4, we can always choose 
i?* to be 

Pi ® . . . ® P„-i ® Qi ® . . . ® Qm-i ®{Ue{Pi®...® P„_i ® Qi ® . . . ® Qm-l)) 

Example 3. Let Ai = [0, a), ^2 = [0, l]\Ai,Bi = [0, 6), P2 = [0, 1] \ ^i, all seen 
as hybrid sets. Let 

[2 < .T < a , , ^ ib Q<x <b 

fix) = < and g[x) = < 

^ ' |0 a<a;< 1 ^ ^ [7 6 < a; < 1 

Wc choose the hybrid (symbolic!) partition Ai, BiQAi, B2, which simultaneously 
refines both. Then after a few computations we get 

/ * 5 = fl2 * 5|}^^ ®* p * Sfr""^®^^ ®* flO * Tfr""^ (2) 

regardless of whether a < 6 or a > 6; in fact either (or both) could be outside of 

[0, 1) and the result, interpreted as a hybrid function, are still correct. Moving 
from one choice of partition to another is done by undoing the distribution, 
performing the change of partition, and using commutativity and associativity 
to regroup like terms. Note that we should have written 2*5 as (a; 1— >■ 2)*{x 1— >■ 5), 
but we chose the above for greater clarity. 

It should be very clear that expressions such as equation 2 are a formal repre- 
sentation of a piecewise function, and need to be interpreted properly in each 
context. 



5.2 Identities for invertible operators 



When the binary operation we use is invertible, we can perform the operations 
at any time, as operands may later be removed from a cumulative result by 
applying their inverses. This may lead to considerable efficiency improvements; 
even though values and their inverses are cancelled in the calculation (leading 
to no net effect) this may be more efficient than retaining an "un-evaluable" 
expression as a pseudo-function. 

Proposition 9. Let be a hybrid function over S where (S*, *) has the struc- 
ture of an Abelian group (where we will use e for the unit and — for the inverse), 
then for all x G supp P, 

{r*f-ni^) = =P(x){|(x,e)i|} 

where —f denotes x i— > —f{x). 

Both equalities follow readily from the definitions. 

5.3 Identities for linear operators 

Definition 20. For a linear operator L, and f^ a hybrid function, 

L{fP) ::= L{x ^ P{x) ■ f{x)) 

Note L{f^) may not be defined even when L{f) is. If the multiplicity function 
P{x) is uniformly bounded, then it will exist. 

Proposition 10. Let f^ be a hybrid function, C[n\Pi « partition of P such that 
each Pi{x) is uniformly bounded, then 

n 

L(n=^L(r') 

i=l 

The above is the fundamental reason why, under Karr's definition, the summa- 
tion identities of the introduction hold. 

Corollary 1. For all total functions f : Z ^ G with G an Abelian group, and 
all £,m,n G Z, 

E = E + E /w- 

i^i<n (^i<m m^i<n 

6 Examples 

We present two examples of the application of hybrid functions to symbolic 

computation problems. The first example is concerned with the arithmetic of 
symbolic matrices, the second presents the idea of merging symbolic spline func- 
tions. 



6.1 Matrix Addition 



Earlier work [10, 11] introduced The idea of support functions has been intro- 
duced previously to represent symbolic matrices — matrices given in terms of 
symbolic regions with undcrspccificd elements and symbolic dimensions and 
defined arithmetic operations between them. The paper [10] presented a sup- 
port function based on half-plane constraints that enables full arithmetic, but 
that suffered from a combinatorial explosion in the number of terms needed to 
express sum or product matrices. The paper [11] moved to a support function 
based on interval addition that automatically dealt with cancellation for neg- 
ative intervals. While this avoided the combinatorial explosion, the approach 
was restricted to certain types of regions and could not be easily generalised to 
matrix multiplication. Hybrid functions and generalised partitions solve both of 
these problems simultaneously. We demonstrate this with the example of matrix 
addition of two 2x2 symbolic block matrices. Let 

where Mi and M2 are nxm matrices, Ai and A2 are of dimensions hi x ki and 
/12 X ^2 respectively, n, m, /ii, /12, fci, ^2 € Nq. 

Let U = I 1 ^ i ^ n A 1 ^ J ^ m} be the set of all cell points in the 

matrices. We define the region occupied by a matrix block similarly, and refer 
both to a matrix block and to the region it occupies by the same name, relying 
on context to distinguish them. We can thus write 

Ml = A^' ® Sf 1 ® Cf ^ ® £)f 1 , M2 = A^^ ® = ® C^^ ® £»f = (3) 

To calculate Mi -|- M2 we: (1) Choose a suitable generalised partition of 
U. (2) Rewrite each block of each matrix into terms restricted to the chosen 
partition. (3) Substitute into the expressions for the matrices. (4) Add the two 
matrices region-wise. 

As established in the section 4, the maximal number of partitions required 
in our case is 4 -|- 4 — 1 = 7. We therefore choose 6 independent regions to 
be Ai,Bi,Ci,A2,B2,C2 and obtain the seventh. Pi, by subtracting all other 
regions from the universe U, 

Pi = [/ e (^1 ® Bi ® Ci ® A2 ® S2 ® C2). (4) 

We can then express regions Di and D2 in terms of Pi: Di = UQ{Ai(BBi(BCi) = 
Pi ® ^2 ® -B2 ® C2, and Da = Pi ® Ai ® Bi ® Ci. 
Rewriting Mi,M2 from Eq. (3), we get: 

Ml = ® Pf 1 ® Cf ^ ® i^fi®^^®s.ec. 

= P)fi ® Af' ® Pf ^ ® Cf ^ ® D^^ ® £>f = ® £)f ^ 

M2 = D^' ® D^' ® D^' ® D^' ® Af^ ® Pf ^ ® C?^ 



This lets us express the sum of the two matrices as the following pseudo- function: 



Ml + M2= (Di + D^f' CD+ (Ai + D^)^' 0)+ [Bi + D^f' ®+ (Ci + D^f^ 

®+ + A2t^ ®+ (Di + ®+ (Ci + Ca)^^ (5) 



Observe that the seven terms of the hybrid function fully capture the result of the 
matrix addition independent of the order of the symbolic dimensions hi,h2, fci, /c2 
of the original blocks. We demonstrate this by evaluating the function for a 
couple of concrete points in the sum matrix. 

First let hi < /i2 and choose a concrete value of a cell where hi < i < h2 
and 1 < j < ki,k2- The point should therefore be in a region composed of 
elements from Bi and A2. Instantiating the multiplicities in equation (5) ver- 
ifies this. Observe that indeed the only regions with multiplicity 1 are Bi = 
I /ii < i < n A 1 s$ J sC fci } and A2 = {(i, j) | 1 < i < /12 A 1 < j < fcs } 
whereas the multiplicities for Ai, B2,Ci,C2 are all 0. Furthermore, we can 
compute the multiplicity for Pi using equation (4). Since the multiplicity of 
the universe U is always 1 every clement is in this partition — we get 
1- (0 + 1 + 0+1 + + 0) = -1. This then yields the computation below, 
which confirms that our element is indeed in the anticipated region (where we 
write region-wise sets {\{Di + 1)2)"^ |} as {Di + .02)"^ to alleviate notation) 



Now assume that the order of the symbolic dimensions for blocks Ai , A2 is 
reversed and we have /i2 < hi. If we now compute the value of a cell with 
with h2 < i < hi and 1 < j < fci,fc2, we get multiplicity for regions Bi,A2, 
but instead multiplicity 1 for Ai,B2. Again the multiplicity for Pi is —1 and 
equation (5) will yield that our cell is in the Ai + B2 region. 

As a final example we compute cell with /ii,/i2 < i < n and ki,k2 < 
j < TO, which is a point from Di and I?2- This time equation (5) simplifies even 
more quickly, as the multiplicities for Ai, A2, Bi, B2,Ci, C2 are all and we only 
have to consider the multiplicity for Pi which is 1 with the same considerations 
as above, yielding Di + D2 as the only term that does not vanish. 

6.2 Symbolic Spline Interpolation 

Numerical computation and manipulation of splines is well understood. Here we 

show how hybrid domain decomposition can be used to support the previously 
unexplored symbolic manipulation of splines. 

Let a = xq < xi < ■ ■ ■ < Xn-i < = 6 be a partition of the interval [a, b] C 
M. We call the Xi knots and assume that for each knot we have a corresponding 
value Ui- Then we can define a spline function S over [a, b] piecewise as 



{Di + D2)"' ®+ {Bi + D2f ®+ {Di + A2f = Bl CD+ Al 



B1+A2 



So{x) X e [xo,xi] 



S{x) = < 



(6) 



^Sn-l{x) X e [Xn-l 




While spline interpolation is traditionally defined for numerical values of the 
Xi,yi pairs, we will now define a symbolic spline function. Let a = Cq < Ci < ■ ■ ■ < 
c„_i < c„ = 6 be symbolic or "abstract" knots with associated symbolic values 
do,di, . . . ,dn, where a di is generally given as a function in Cj. We define spline 
segments Sci,ci+i (x) for i = 0, . . . , n — 1. That is, the segments are parameterised 
with respect to two knots and their values. Let P = Pi (B ■ ■ ■ (B Pn with Pj = 
[ci-i,Ci] be a generalised partition. We then define a symbolic spline function 
S{x) = S^^^^{x) ® • • • ® Sc^_^^cn{^)^"' ■ For clarity we often omit the (x) part of 
the term. 

Define the merge of two spline segments n S^y to be 5rSa,a'),mm(fc,fe'r 
Clearly this merge will be empty if the two segments do not overlap, otherwise it 
will be the smallest possible spline for the overlapping interval of the segments. 

Now let Pi © • • • © P,i and Qi ® • • • ® Qm be two generalised partitions of the 
interval [a, 6] and let S = S^^\^^(D- ■ ■(DS^^_^^^^ and T = T^^^^^ ©• • -©T^^"^^ ,^^ be 
two symbolic splines. Observe that the di here are knots and not knot values. We 
can then define the merge 5 N T as a binary operation on two hybrid functions 
as given above in proposition 8. 

We observe the merge operation using a simple example!. Lcrt P = PiffiP2 and 
Q = Qi®Q2 be the generalised partitions a < c < h and a < d < b of our universe 
[a, b], respectively. Let S = S^^^CDS^^,^ and T = T^}(DT2^ be two symbolic splines. 
We choose a common refinement as Pi,Qi,R — U Q (Pi © Qi). We can then 
write S = S^l, © S^^"^' = S^^, © 5« © Sf^l and similarly T = T^;, © T^^ © r«,. 
When we merge both symbolic splines we get 

5 X T = {Sa,c N Ta,bf' ©« (5e,6 M T„,d)«i ©^ {Sc,b X (7) 

If we now fix the order of our symbolic knots to be a < c < d < 6 we can 
evaluate the three components of our spline. First let a < x < c, which means 
-Pi = Qi = 1 and P = — 1 and (7) evaluates to Sa,c x Ia,d which yields a spline 
between knots a, c. Similarly the other two segments evaluate Sc,b ^ T^^d and 
Sc,b ^ ^d,6) which yields splines for the intervals [c, rf] and [d,h\, respectively. 

7 Conclusion 

We have presented a framework of generalised partitions and domain decompo- 
sition based on hybrid sets. This framework has a number of pleasing properties 
and unifies a number of ad hoc notions in common use. More importantly for our 
purposes, this representation allows easy manipulation of and reasoning about 
partitions whose pieces are defined symbolically. These specialise correctly for 
all choices of parameters by negative and positive multiplicities cancelling as 
needed. The representation (see for example equation (2)) needs to be carefully 
interpreted, as out of context simplification can result in meaningless results. 
But if the rules we lay out are followed, our compact representation effectively 
allows one to compute with very general piecewise-dcfined functions. 

Although a number of previous authors have studied hybrid sets, our study 
of functions over hybrid sets, generalised partitions, the superposition © and 



marked superposition (D* operators appear to all be new. Our applications to 
computation and reasoning are certainly new. 

There remain several intriguing directions for future work. These include 
normal forms for piecewise functions defined on hybrid partitions, simplifica- 
tion of intermediate expressions involving linear operators or inverse elements 
and creating partition schemes from algebraic specialisation properties, e.g. in 
Cylindrical Algebraic Decomposition. We have implemented a prototype Maple 
package for hybrid sets and hybrid functions, but a more general implementation 
would be of interest. 
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